cpi uiaii 



Estimation of the effect of hyperonic three-body forces on the 
maximum mass of neutron stars 
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PACS 13 . 75 . Ev - Hyperon-nucleon interactions 
PACS 26 . 60 . Kp - Equations of state of neutron-star matter 
PACS 26.60.-c - Nuclear matter aspects of neutron stars 
PACS 97 . 60 . Jd - Neutron stars 

Abstract. - A model based on a microscopic Brueckner-Hartree-Fock approach of hyperonic 
matter supplemented with additional simple phenomenological density-dependent contact terms 
is employed to estimate the effect of hyperonic three-body forces on the maximum mass of neutron 
stars. Our results show that although hyperonic three-body forces can reconcile the maximum 
mass of hyperonic stars with the current limit of 1.4 — 1.5Mq, they are unable to provide the 
repulsion needed to make the maximum mass compatible with the observation of massive neutron 
stars, such as the recent measurements of the unusually high masses of the millisecond pulsars 
PSR J1614-2230 (1.97 ± O.O4M ) and PSR J1903+0327 (1.667 ± 0.021M©). 



Neutron stars are the remnants of the gravitational col- 
lapse of massive stars during a Type II, lb or Ic supernova 
explosion. Their masses and radii are typically of the or- 
der of 1 - 2M Q (Mq ~ 2 x 10 33 g is the mass of the Sun) 
and 10 — 12 km, respectively. With central densities in 
the range of 4 — 8 times the normal nuclear matter satu- 
ration density, eo ~ 2.7 x 10 14 g/cm 3 (p ~ 0.16 fm~ 3 ), 
neutron stars are most likely among the densest objects 
in the universe pQ. Nowadays, it is still an open question 
which is the true nature of neutron stars. Traditionally 
the core of neutron stars has been modeled as an uniform 
fluid of neutron rich nuclear matter in equilibrium with 
respect to the weak interaction (/3-stable nuclear matter) . 
Nevertheless, due to the large value of the density, new 
hadronic degrees of freedom are expected to appear in ad- 
dition to nuclcons. Hyperons, baryons with a strangeness 
content, are an example of these degrees of freedom. Since 
the pioneer work of Ambartsumyan and Saakyan [2] the 
presence of hyperons in neutron stars has been studied by 
many authors using either phenomenological [3HS] or mi- 
croscopic [THlOj approaches. Hyperons may appear in the 
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inner core of neutron stars at densities of about 2 — 3po- 
At such densities, the nucleon chemical potential is large 
enough to make the conversion of nucleons into hyperons 
energetically favorable. This conversion relieves the Fermi 
pressure exerted by the baryons, and makes the equation 
of state (EoS) softer. As a consequence, the maximum 
mass of the star is substantially reduced to values that, 
in microscopic calculations [Tl-I10j. can be even below the 
"canonical" one of 1.4 — 1.5M Q inferred from precise neu- 
tron star mass determinations [TT]. It has been suggested, 
see e.g. Ref. [T2], that a possible solution to this problem 
is the inclusion of three-body forces (TBF) involving one 
or more hyperons (i.e., nucleon- nucleon- hypcron (NNY), 
nucleon- hyperon-hypcron (NYY), and hyperon-hyperon- 
hypcron (YYY)). These forces could eventually provide 
the additional repulsion needed to make the EoS stiffcr 
and, therefore, the maximum mass compatible with the 
current observational limits. In this letter, we want to 
establish numerical lower and upper limits to the effect 
of these forces on the maximum mass of neutron stars 
that can guide more sophisticated calculations in the fu- 
ture. To such end, we use a model based on a micro- 
scopic Brueckner-Hartree-Fock (BHF) approach of hyper- 
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onic matter [7T-|9"I[T5] supplemented with additional simple 
phenomenological density-dependent contact terms that 
account for the effect of nucleonic and hyperonic TBF. 

Our calculation starts by constructing all the baryon- 
baryon G matrices, which describe the interaction between 
two baryons in the presence of a surrounding medium. The 
G matrices are obtained by solving the coupled-channel 
Bethe-Goldstone equation, written schematically as 



G(ui)B 1 B 2: B 3 B i — Vb 1 B 2 ,B 3 B 4 + ^ Vb 1 B 2 ,B z B ] 
QBiBj 



B,B, 



w - E Bi - E Bj + it] 



G{w) Bi B j ,B 3 B i , (1) 



where the first (last) two subindices indicate the initial 
(final) two-baryon states compatible with a given value S 
of the strangeness, namely nucleon-nucleon (NN) for S = 
and hyperon-nucleon (YN) for S = — 1; V Bi b 2 ,b 3 b 4 is 
the bare baryon-baryon interaction (NN— s-NN, AN— »AN, 
AN->£N or £N->£N ); Q BiBj is the Pauli operator, that 
prevents the intermediate baryons Bi and Bj from being 
scattered to states below their respective Fermi momenta; 
and the starting energy u> corresponds to the sum of the 
nonrclativistic single-particle energies of the interacting 
baryons. 

The single-particle energy of a baryon Bi is given by 



E Bi {k)=M Bl 



2 1,2 



2M r 



Re[U Bl {k)} 



(2) 



Here denotes the rest mass of the baryon, and the 
single-particle potential Ust represents the average field 
"felt" by the baryon owing to its interaction with the other 
baryons of the medium. In the BHF approximation, U Bi 
is calculated through the "on-shell energy" G-matrix, and 
is given by 



3 J k' 



x(kk'\G(E Bi {k) + E Bj (k , )) Bt B i ,B i B i \kk') A , (3) 

where n Bj (\k\) is the occupation number of the species 
Bj, and the index A indicates that the matrix elements 
are properly antisymmetrized when baryons Bi and Bj 
belong to the same isomultiplet. We note here that 
the so-called continuous prescription has been adopted 
for the single-particle potentials when solving the Bethe- 
Goldstone equation, since, as shown by the authors of Ref. 
P3] j the contribution to the energy per particle from three- 
hole line diagrams is minimized in this prescription. All 
the calculations carried out in this letter have been per- 
formed with the realistic Argonne V18 [15] NN force and 
the Nijmegen soft-core NSC89 [16] YN interaction. We 
recall that the NSC89 potential does not contain hyperon- 
hyperon (YY) components. 

Once a self-consistent solution of Eqs. (H])-© is ob- 
tained, the energy density can be calculated in the BHF 



approximation as 
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TBF can be introduced in the BHF approach by adding 
effective density-dependent two-body forces to the baryon- 
baryon interactions V when solving the Bethe-Goldstone 
equation for the G-matrices. These effective forces are ob- 
tained by averaging genuine baryon-baryon-baryon TBF, 
Wz{ri,fj,rk), over the third baryon coordinates [17], 



(5) 



where Tij is the relative coordinate of baryons Bi and Bj , 
and n[fi , fj , rj.) is an appropriate three-body correlation 
function. This procedure has been the usual one to in- 
troduce nucleonic TBF in BHF and other non-relativistic 
many-body approaches of nuclear matter (see e.g., Ref. 
[Tg]). Nevertheless, since there is an almost complete lack 
of experimental data, and few theoretical studies [T9"H2"T] . 
on the genuine hyperonic TBF, instead of using the aver- 
age procedure described, in this letter we adopt an alter- 
native strategy. First, we construct, as described above, 
the hyperonic matter EoS within the BHF approach us- 
ing only two-body NN and YN forces, and then, we add 
to it simple phenomenological density-dependent contact 
terms that account for the effect of both nucleonic and 
hyperonic TBF. The reader could think that in doing that 
the predictive power of our model is partially lost. How- 
ever, we would like to stress once more that the aim of this 
work is to perform a simple estimation of the effect of hy- 
peronic TBF on the maximum mass of neutron stars that 
can guide and motivate more realistic and sophisticated 
studies of such forces and their effects. We want also to 
point out that by using a simple parametrization of the 
hyperonic TBF we can easily explore whether or not such 
forces can make the EoS stiff enough so that the maximum 
mass of hyperonic stars can be reconciled with the obser- 
vational limits. Following Balberg and Gal (see Eqs. (2) 
and (3) of Ref. [6]), we assume these density-dependent 
contact terms to have a Skyrme-like form 
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where x and y denote any baryon species. Note that, to 
be precise only the second term actually accounts for the 
effect of three- and multi-body forces which are responsi- 
ble for the repulsion that dominates at high densities. The 
first one is in fact a two-body term that modifies slightly 
the BHF part of the energy density by yielding additional 
attraction at low densities. It can be considered as a rear- 
rangement of the two-body contribution due to the inclu- 
sion of TBF. It mimics the modification that the effective 
density-dependent two-body force, V B B , would introduce 
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on the G-matrices. Then, we have the following contact 
terms contributing to the energy density 



_7Niv + l 



o>nnPn + vnnPn 

a ANPA p N + b AN p APN ( ^ + Pn 
\ PA + Pn 

+ a^ N p^p N + b^NP^PN 
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fig* 
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where p N = p n + Pp , p s = p s - + ps a + Ps+, and we 
assume charge independence. We note that the first two 
terms mime the contribution of the effective two-body in- 
teration, VJyjy, obtained when the forces NNN and NNY 
are averaged over the coordinates of one of the nucleons 
and the hyperon, respectively (see Eq. (0). Similarly, 
the remaining terms mimic the contributions of the effec- 
tive Vy J forces built from the average of the NNY and 
NYY ones over the coordinates of the appropriate nucle- 
ons and hyperons. Note that the NYY forces can also 
contribute, in addition to the YYY ones, to the effective 
density-dependent YY forces if the average is taken over 
the nucleon coordinates. For consistency with the use of 
the NSC89 YN potential, which as it is said, does not de- 
scribe any YY interaction, none of the terms of Eq. (|7|) 
mimic neither this other contribution of NYY, nor the 
ones of YYY. The complete baryonic energy density is 
then simply obtained by adding cct to the BHF energy 
density of Eq. Q. 

The parameters a^N, and jnn are fitted to repro- 
duce the density, po, energy per particle, E/A, and incom- 
pressibility coefficient, K 00l of symmetric nuclear matter 
at saturation. We choose the saturation point to be at 
po = 0.16 far 3 and E/A = -16 MeV. Due to the uncer- 
tainty still existing on the precise value of i^oo, we con- 
sider four values of jnn between 2 and 3.5, giving 
in the range 211 — 285 MeV, compatible with the value of 
Koo = 240 ±40 MeV supported nowadays by experimental 
data [28]. The values of cinn, bNN, 7 aw and -ftToo are given 
in Table [TJ We note that due to the charge independence 
assumed in Eq. (JJJ neither the symmetry energy, E sym (p) , 
nor its derivative parameter L = 3podE sym (p) / dp\ Po will 
be affected by the inclusion of the contact terms. Their 
values at p n = 0.16 fm~ 3 are 28.3 and 46.6 MeV, respec- 
tively. 



Inn 


a-NN 


bNN 






[MeV fm 3 ] 


[MeV firr 7 ^] 


[MeV] 


2 


-33.44 


213.02 


211 


2.5 


-22.08 


355.03 


236 


3 


-16.40 


665.68 


260 


3.5 


-12.99 


1331.36 


285 



Table 1: Values of the parameters ajvjv, 6jvjv, 7jvjv, and the 
incompressibility coefficient K^. 



TBF. There is no physical reason for which these param- 
eters should scale with the same factor x with respect to 
their nucleonic counterparts. However, we have checked 
that allowing different scaling factors for clyn and by^, 
we always obtain maximum masses within the range of 
values shown in Table [5J As an example, we note that for 
7AfAr=3.5 we obtain M max = 1.61M for ayN = o<nn/3 
and by n = 6 aw, which is very similar to the value of 
I.6OM0 shown in Table [2] for x = 1. In view of this, 
and for simplicity, we keep the same scaling factor x for 
both ay 7v and byN- Finally, the value of —28 MeV for 
the binding energy of a single A in nuclear matter at sat- 
uration, extracted from the extrapolation of hypcrnuclear 
data [29], is then used to determine the parameter ^yn 
through the relation 



= -28 MeV = U A (k = 0) + a YN pa + b YN pl 



7YN 



(8) 

where U A (k = 0) (—30.8 MeV in our calculation) is the 
BHF single-particle potential of the A (see Eq. ([3])) at zero 
momentum, and the last two terms account for the contact 
term contribution. 

Once the total energy density e = (bhf + £ct + (l 
{cl being the contribution of noninteracting leptons), is 
known, the composition and the EoS of neutron star mat- 
ter can be obtained from the requirement of equilibrium 
under weak interaction processes, pn = bip n — qip e (pi and 
qi denoting the baryon number and charge of species i) 
and electric charge neutrality, J^iQiPi = 0- The chem- 
ical potentials of the various species and the pressure 
are computed from the usual thcrmodynamical relations, 
Pi = de/dpi and P = p 2 d(e/ p) /dp. Finally, knowing 
the EoS, the equilibrium configurations of static neutron 
stars are obtained by solving the well-known Tolman- 
Oppcnheimcr-Volkoff (TOV) equations [T|. 



In order to reduce the number of parameters, in this ex- 
ploratory work, we assume for simplicity that TBF involv- 
ing A and E hyperons are the same. Therefore, we take 
a A N = &T.N-, b A N — bsN and 7a/v = Jsn- The parameters 
fly n and byN (Y = A, E) are considered free. We take 
them to be a fraction x of the parameters ajy/v and 6atat 
(ayN = xclnn and byN = x bNN with x = 0, 1/3, 2/3 and 
1) in order to explore different strengths of the hyperonic 



Table [5] shows the maximum neutron star mass, central 
baryon number density and central speed of sound for dif- 
ferent values of the contact term parameters. The range 
of variation of these parameters allows us to explore dif- 
ferent EoS and establish a lower and an upper limit to 
the effect of hyperonic TBF on the maximum mass. One 
can get an idea of the relative importance of these forces 
on the EoS (and, therefore, in some way also of the rel- 
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Inn 


X 


1YN 




Pc 


V s 







- 


1.27(2.22) 


1.35(1.07) 


0.46(1.03) 




1/3 


1.49 


1.33 


1.33 


0.48 


2 


2/3 


1.69 


1.38 


1.29 


0.52 




1 


1.77 


1.41 


1.24 


0.54 







- 


1.29 (2.46) 


1.19(0.92) 


0.43(1.17) 




1/3 


1.84 


1.38 


1.16 


0.49 


2.5 


2/3 


2.08 


1.44 


1.12 


0.54 




1 


2.19 


1.48 


1.09 


0.56 







- 


1.34(2.72) 


0.98 (0.79) 


0.40(1.34) 




1/3 


2.23 


1.45 


0.97 


0.50 


3 


2/3 


2.49 


1.50 


0.94 


0.55 




1 


2.62 


1.54 


0.90 


0.58 









1.38 (2.97) 


0.87 (0.69) 


0.38(1.47) 




1/3 


2.63 


1.51 


0.86 


0.51 


3.5 


2/3 


2.91 


1.56 


0.83 


0.56 




1 


3.05 


1.60 


0.80 


0.59 



Table 2: Maximum neutron star mass, central baryon number 
density and central speed of sound for different values of the 
contact term parameters. The results for x = correspond to 
the case when only nucleonic TBF are considered. In brakets 
are given the corresponding values when the presence of hyper- 
ons is neglected. Masses are given in Mq whereas the central 
baryon density, p c , is given in fm~ 3 and the central speed of 
sound is given in units of c. 
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ative weight of the phenomcnological part of our calcu- 
lation with respect to the microscopic one), by evaluat- 
ing the ratio between the last four terms of Eq. ([7]) , and 
the sum Es I Efc"s.(|fc|)t / s I (fc)/(2V) (see Eq. ©). We 
find that, in average, this ratio is smaller than 0.2 for 
p < 4po, it ranges between 0.2 and 0.5 for densities up to 
~ 5/?o, and it is larger than 1 for p > 6pg, clearly showing 
that the relative importance of hyperonic TBF increases 
for larger densities. The results for x = correspond to 
the case when only nucleonic TBF are considered {i.e., 
clyn = byN — 0). In brakets arc given the correspond- 
ing values when the presence of hyperons is neglected in 
the EoS. Note that in this case the resulting maximum 
mass is relatively large, ranging from 2.22M Q for ^nn = 2 
to 2.97M Q for 7wjv = 3.5. The presence of hyperons in- 
duces here a reduction of the mass to values in the interval 
1.27 - 1.38M Q , well below the limit of 1.4 - 1.5M . Note 
that the range of variation of the maximum mass is about 
O.IIMq in this case, compared to a range of ~ 0.75M Q 
when hyperons are absent. This is a consequence, as al- 
ready pointed out in Rcf. [S], of a strong compensation 
mechanism caused by the appearance of hyperons which 
makes the maximum mass quite insensitive to the pure nu- 
cleonic part of the EoS. As expected, the central density 
decreases when increasing the effect of three-body forces 
(the pressure is larger and consequently the object is less 
compact), but at the same time the speed of sound in- 
creases, because the EoS is stiffcr. We note that when 
the presence of hyperons is neglected the EoS is always 



Fig. 1: (Color online) Upper panel: /3-stable matter EoS. Lower 
panel: Mass-radius relation for different EoS. Circles indicate 
the central baryon number density, central pressure, mass and 
radius of the maximum mass stellar configuration. Horizontal 
lines show the masses of the pulsars PSR J1614-2230 [3D], PSR 
J1903+0327 [31] and the Hulse-Taylor one [32]. See the text 
for details. 



supraluminical. This is not surprising, since our approach 
is a non-relativistic one, and causality is, therefore, not 
guaranteed. However, note that as soon as hyperons are 
present in matter, the softening of the EoS induced by 
their presence is such that in these cases the EoS remains 
always causal. It is clear from Table[5]that hyperonic TBF 
provide additional repulsion making the EoS stiffer, and 
the maximum mass larger. For a fixed value of the ex- 
ponent jnn the maximum mass increases when increas- 
ing x (i.e., ~/yn)- This is an expected result, since by 
increasing x we are increasing the strength of the hyper- 
onic TBF and, as a consequence, the EoS becomes stiffcr, 
and the maximum mass larger. We have checked that 
the rate of increase of M max with x is slightly quadratic. 
The stiffer EoS including hyperonic TBF is obtained for 
Inn = 3.5 and x = 1 (jyn = 3.05), and allows for a max- 
imum mass of about 1.60M©. We note that although the 
inclusion of hyperonic TBF can reconcile the maximum 
mass of hyperonic stars with the "canonical" value, they 
are, however, unable to make the maximum mass compat- 
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ible with the observation of massive neutron stars, such as 
the recent measurements of the unusually high masses of 
the millisecond pulsars PSR J1614-2230 (1.97 ± 0.04M Q ) 
[30] and PSR J1903+0327 (1.667 ± O.O21M ) [31]. One 
could in principle increase arbitrarily the strength of the 
hypcronic TBF by increasing the value of the fraction x 
in order to get a very stiff EoS which supports such large 
masses. However, in our opinion values of x > 1 give rise 
to EoS unrealistically stiff, the reason being the following: 
we know that the strength of the two-body YN interaction 
is smaller than that of the NN one (e.g., the single-particle 
potential of a A in symmetric nuclear matter at saturation 
for zero momentum is about 1/3 that of the nucleons). 
Therefore, it is quite natural to think that probably the 
strength of hyperonic TBF is either smaller or as large 
as the pure nucleonic ones, but not larger. Although this 
statement seems a bit speculative, and only a more realis- 
tic and sophisticated determination of hyperonic TBF can 
give a definite answer, we consider the value of 1.60M© 
a reasonable upper limit for the maximum mass of neu- 
tron stars with hypcronic TBF. This upper limit could be 
slightly larger, not more than 5%, if we would have taken 
a^N ^ a an, few ^ b\N, and 7ew ^ Jan, and assumed a 
slightly repulsive single-particle potential for the E~ in nu- 
clear matter at saturation density. We point out, however, 
that the uncertainties related with the binding-energy of 
the £~ in nuclear matter at saturation density are still 
very large (see e.g., Refs. [35H36"] ). Clearly, 1.27M Q is the 
lower limit corresponding to the softer EoS. 

Fig.rjjshows as a summary the EoS (upper panel) for the 
stiffcr and the softer pure nucleonic (curves 1 and 2) and 
hypcronic (curves 3 and 4) stars, and their corresponding 
mass-radius relation (lower panel). Curves 1 and 2 show 
respectively the results for jnn = 3.5 and jnn = 2 with- 
out hyperons, while curves 3 and 4 correspond to the cases 
Jnn = 3.5, x = 1, and jnn = 2,x = 0. Circles indicate 
the value of the central baryon number density, central 
pressure, mass and radius of the maximum mass stellar 
configuration. The horizontal lines show the masses of the 
pulsars PSR J1614-2230, PSR J1903+0327 and the Hulse- 
Taylor one (1.4414 ± 0.0002) [32]. The strong softening of 
the EoS due to the presence of hyperons and the conse- 
quent reduction of the maximum mass is clearly seen. The 
maximum masses of hyperonic stars lay in a narrow range 
from 1.27 to 1.60M©, that is still compatible with the 
mass of Hulse- Taylor pulsar, but is well below the masses 
of PSR J1903+0327 and PSR J1614-2230. 

Summarizing, we use a model based on a microscopic 
BHF approach of hypcronic matter supplemented with ad- 
ditional simple phenomcnological density-dependent con- 
tact terms to establish numerical lower and upper limits 
to the effect of hyperonic TBF on the maximum mass of 
neutron stars. Assuming that the strength of these forces 
is either smaller or as large as the pure nucleonic ones, our 
results show that maximum masses of hyperonic stars lay 
in a narrow range from 1.27 to 1.60M Q which is still com- 
patible with the "canonical" value of 1.4 — 1.5M Q , but 



it is incompatible with the observation of massive neu- 
tron stars, such as the recent measurements of a mass of 
1.97±0.04M Q for the millisecond pulsar PSR J1614-2230, 
and a mass of 1.667 ± 0.021A/ Q for the PSR J1903+0327 
one. We hope that this exploratory work can serve as 
a motivation to perform more realistic and sophisticated 
studies of hyperonic TBF and their effects on the neutron 
star structure, since they have the last word on this issue. 
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